Entropic phase separation of linked beads 
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We study theoretically a model system of a transient network of microemulsion droplets connected 
by telechelic polymers and explain recent experimental findings. Despite the absence of any specific 
interactions between either the droplets or polymer chains, we predict that as the number of poly- 
mers per drop is increased, the system undergoes a first order phase separation into a dense, highly 
connected phase, in equilibrium with dilute droplets, decorated by polymer loops. The phase transi- 
tion is purely entropic and is driven by the interplay between the translational entropy of the drops 
and the configurational entropy of the polymer connections between them. Because it is dominated 
by entropic effects, the phase separation mechanism of the system is extremely robust and does not 
depend on the particlular physical realization of the network. The discussed model applies as well 
to other systems of polymer linked particle aggregates, such as nano-particles connected with DNA 
linkers. 
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PACS numbers: 

Introduction: Equilibrium or transient gels are 
network-forming systems, examples of which include sur- 
factant solutions, gels of biological molecules or synthetic 
polymers [J II SB II IS II GB 111 III Understanding 
the phase behavior and structure of such systems is an 
active field of research 0, E Q E E El C1M III , 
and has a wide range of practical applications [lLI3.l7lllfjj. 
A particularly elegant experimental realization of a tran- 
sient network has been reported in |2l| . The system con- 
sists of oil-in- water microemulsion droplets (which we call 
either drops or beads) connected by telechelic polymers 
(see Fig. |T]b)); the latter have a hydrophilic backbone 
with a hydrophobic group at each chain end. Mixtures 
of telechelic polymers and emulsions have a wide range 
of technological applications, including paints, cosmetics 
and enhanced oil recovery. Precision control of the struc- 
tural and rheological properties of materials is essential 
for achieving good performance. Such control is currently 
achieved using the telechelic additives, that form a tran- 
sient network with controlled rheological and structural 
properties |2l| . Apart from its high applicative interest, 
the telechelic-microemulsion mixtures serve as a model 
system for a general class of transient networks. The 
advantage of this particular system is that the parame- 
ters that control the thermodynamics and structure can 
be easily identified and independently controlled: the 
concentration of possible nodes (the droplets) and the 
connectivity of the network (the number of polymers per 
droplet). This is in contrast, for example, with binary 
mixtures of telechelics, where one cannot separately con- 
trol the number of nodes, formed by the associating chain 
ends. A similar system of gold nano-particles connected 
by short DNA segments, has been of a considerable in- 
terest recently HQ 

Due to the high hydrophobic energy of the stickers (~ 
10 — 20kftT ), the chain ends are constrained to lie within 
the droplets, and the number of dangling ends outside the 
droplets is statistically insignificant. Nevertheless, the 



polymer ends can detach from a droplet and switch be- 
tween loop (with both ends inside the same droplet) and 
bridge (with the chain ends residing in different droplets) 
configurations when the droplets are close enough. Poly- 
mer ends can also be exchanged between the droplets 
during droplet collisions. Rheological experiments show 
that the stress relaxation times, which are related to the 
timescale for chain rearrangement between the droplets, 
are of the order of one second [2lj. Thus, the timescales 
of the observation (days) are much larger than the time 
for the polymers to explore different connectivity config- 
urations, and the system can be considered to be in equi- 
librium. Although the system is athermal, it has been 
found J21j that for high polymer to bead ratios, the sys- 
tem undergoes a first-order phase separation. Indepen- 
dently, a non-thermodynamic, structural transition was 
observed, where an infinite connected network is formed, 
as reflected in rheological measurements. The experi- 
mentally observed phase diagram is shown in Fig. ^a). 

In this Letter, we investigate theoretically the phase 
behavior and structure of this experimental system. We 
predict that as the average number of polymers per bead 
is increased, the system undergoes a first-order phase sep- 
aration into a highly connected dense network that coex- 
ists with a dilute phase of disconnected beads, decorated 
with polymer loops. This explains the experimental re- 
sults of Ref. [2l| and predicts how the polymer properties 
control the phase behavior. The predicted phase separa- 
tion has a purely entropic origin; there are no energetic 
interactions among the polymers or droplets. The phase 
separation occurs because the loss of the translational 
entropy of the droplets is overcompensated by the high 
configurational entropy of the polymer connections in the 
dense network. The analytical results are confirmed by 
Monte Carlo simulations that also predict the phase sep- 
aration. 

Analytical model: We first discuss the thermodynamic 
behavior of the system. A total of N beads and N p 
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Z\ = qiNe~ e + cN J2r 1r e~ En ^ T where e is a free energy 
cost (in units of fceT) of a looped polymer, that measures 
the entropic cost of both ends being confined to the same 
droplet. Taking into account the indistinguishability of 
the polymers, the configurational partition function of 
the total N p polymers, that determines the polymer dis- 
tribution among the beads is thus 
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FIG. 1: a)Description of the connected microemulsion. The 
telechelic polymers can either link two oil droplets or loop on 
a single one. b) Experimentally observed phase behavior of 
the connected microemulsion as a function of the bead density 
c and the mean number of stickers per droplet, r. 



polymers are distributed in space so that a polymer ei- 
ther connects two beads or loops on a single bead. The 
total free energy of the system, which in this athermal 
system amounts to the configurational entropy, includes 
two contributions. The first, is the translational en- 
tropy of the beads which we take as the lattice-gas en- 
tropy of mixing of a dispersion of hard spheres, Sq(c) — 
— (c In c+ (1 — c) In ( 1 — c) ) , where c is the volume fraction 
of the beads and So(c) is the entropy per site on the lat- 
tice. One should bear in mind that c is the effective vol- 
ume fraction of the droplets, together with the thetered 
polymers. The experimentally controlled density of pure 
R 3 

beads is roughly C (R G+ °R yi ■ The theory treating the ef- 
fects of the excluded volume of the polymers, will be 
presented elsewhere. The second contribution to the free 
energy of the system is the configurational entropy of 
distributing the polymers among the beads. A polymer 
can either loop onto a single bead or bridge between two 
droplets. For a single polymer, there are qiN available 
looped states, where N is the number of the drops in 
the system and qi is the number of positions available to 

a sticker of size a at the surface of a drop: qi ~ 4 J^|° . 
A bridging polymer connecting two drops at a distance 
R, is stretched to a length R, with an entropic cost (in 
units of k-g,T and assuming ideal chains, for simplicity) 
Er ~ |i? 2 /i?g where Rq is the polymer radius of gy- 
ration. From simple geometric considerations, the aver- 
age number of droplet pairs separated by a distance R 
is hcN ^§ g(R) = qRcN where Rq is the radius of a 
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drop and g (R) is the normalized density-density correla- 
tion function of the drops. As a first step, in the mean 
field approximation we take g(R) = I for R > Rq. As 
we show later, the average number of looped polymers 
per bead is small. Thus, because there is no restric- 
tion on the number of polymers connecting two given 
drops (except at very high polymer concentrations, out- 
side the range of the experiments), the polymers are in- 
dependent. The partition function of a single polymer is 



where q = f dR 47T Q e Er / t , is the effective number of 

droplet pairs. It is determined by the balance of the 
Boltzmann factor for the stretching energy, that strongly 
decreases with distance, and the number of available po- 
sitions, 4v S> that increases with distance R. 
In the thermodynamic limit, where the number of drops, 
polymers and the number of lattice sites, V, all become 
infinite, the free energy per site ( in units of k^T), 
f = -ylnZ p - S (c) is 

/ = clnc+(l - c) In (I - c)+(f>(\n(f> - l)-01n (qc 2 + cqieT 

(2) 

where c = N/V and 4> — N p /V are the densities of 
the droplets and of the chains respectively. The first 
two terms account for the translational entropy of the 
droplets and the third term is the translational entropy 
of the polymer chains (the excluded volume of the poly- 
mers is neglected at this stage). The last term is the ef- 
fective, entropically-induced attraction between the poly- 
mers and droplets; this term causes the free energy to 
decrease when the concentration of the beads or poly- 
mers/beads ratio is increased. We stress that there are 
no specific interactions between the droplets and/or the 
polymers, and the resulting thermodynamic attraction 
is of purely entropic origin. The system is thermody- 
namically stable if the free energy / (4>, c) (eq. ©is a 
convex function of both (f> and c. Formally, this means 
that the second variation of the free energy, S 2 F = 
F^5(f) 2 + F cc 5c 2 + 2F c< pS(j)5c is positive so that the sta- 
bility matrix S = ^jj"^ is positive definite i.e. 

possessing two positive eigenvalues. Because F^ = ^ is 

always positive, it is sufficient for the determinant Det S 
to be positive, in order to guarantee convexity. Using Eq. 
© for / (0, c), one finds Det § = -2 ^^ + 

and the system is thermodynamically stable only if 



2(f>/c < 



c + ^-e c 
c(l - c) 



(3) 



where 2<f)/c = r is the average number of stickers (poly- 
mer chain ends) per droplet. This stability condition de- 
fines the spinodal line in the (c, r) plane as shown in Fig. 
121 Note that because the destabilizing effect comes from 
the cross-terms, F c( k, the tie-lines are not horizontal, but 
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connect points with different values of r. If r > c (i-c) ' 
the system is thermodynamically unstable and phase sep- 
arates into a system of dense droplets that are highly con- 
nected by polymers, that coexists with a dilute system 
of almost disconnected droplets, decorated with polymer 
loops. The latter observation stems from the fact that 
the average fraction of the looped polymers A, is given 
by: 



dlii(Zi) 
de 



(4) 



as follows from eq.QJ. Thus, the fraction of looped 
polymers A tends to unity in the dilute phase, where 
From the equality of the polymer chemi- 

in the coexisting phases, 



qc 2 -\-qie~ € c 



cal potential, — In 
it follows that the phases that coexist lie along the lines 
r = m(c + ^-e~ E ) in the (c, r) plane; where m is a con- 
stant determined by /i^, m = e M */g. For any given 
m > 4, this line intersects the spinodal (eq. |3J) at two 
points, as shown in Fig. 2. For m — 4 there is only 

one solution \ (l + ^ e j\ anc ^ *he une ^ s ^ an " 
gent to the spinodal; this determines the location of the 
critical point which is in our mean field model is given 
by c = 0.5 independent of the value of the control pa- 
rameter ■ 



and the critical polymer to bead ratio, 



= 2 



For m < 4 there is no phase coex- 



istence. The predicted phase diagram of the system is 
shown in Fig. 2. Comparing the phase diagram of Fig. 
2 with the experimentally observed one of Ref. [2lJ 

We now estimate the parameter e that reflects the re- 
duction in the entropy of a polymer due to the constraint 
that both ends reside in the same droplet. In a simple 
approximation, the number of configurations available to 
a polymer with radius Rg, with both ends constrained 
to a volume v, that is small relative to the total volume 

R?l 

available, is proportional to -55- ~ -S- where I is the 

R? I 

length of the hydrophobic sticker. Therefore e~ e ~ 

G 

for i? < Rg an d saturates to unity for R > Rg- The 
parameters qi , q and e (which is entropic in origin) , com- 
bine to yield the single control parameter of the system: 
qie~ e I q. Using typical values of the droplet radius, poly- 
mers radius of gyration and the size of the hydrophobic 
stickers from Ref. [2lJ , we estimate the control parameter 
qie- f /q ~ 2 - 3. 

It is important to realize that the predicted phase sep- 
aration cannot be attributed to the direct attraction 
between the beads induced by the connecting poly- 
mers; the total polymer stretching energy site is E ~ 



x J E R e- E R' T q R dR 



R' T q R dR 



and is irrelevant for the phase sepa- 



ration because it is linear in polymer density <j). 

Due to its entropic nature, the phase separation is 
extremely robust and is independent of the detailed as- 
sumptions about the polymer properties. It persists even 
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FIG. 2: Phase diagram of drops connected with polymers. 
The thick line is the spinodal line of the phase separation for 
qie~ e /q = 2. Above this line the system becomes thermo- 
dynamically unstable. The critical point is at c = 0.5 and is 
shown as a black dot. Note that the critical point is not at the 
minimum of the spinodal. The tie-lines are shown as dotted 
lines in the phase separation region. Note that they are not 
horizontal. The dashed line shows the percolation threshold 
calculated for an f.c.c. lattice with q — 16. Below the perco- 
lation line, the system is in the fluid state, while above it a 
connected gel is formed. Onr should remember that the ex- 
perimentally controlled volume fraction of the beads is equal 



to c 



(fl G + iio) a 



in the limit of very short and rigid chains, that connect 
only nearest neighbor drops, complementary to the case 
of flexible chains discussed above. Neglecting the single 
bead loops (which can be shown to be unimportant in this 
limit), the total number of ways to distribute N p indistin- 
guishable polymers among P = ^qcN nearest neighbor 



. • (Ar p +P-1)! 

pairs is A r p! (p_ 1 ) ! 



27| and the total free energy per unit 

volume is / ~ —So (c) + \qc 2 (ln(igc 2 ) — l) — \qc 2 \n.(j> 
when r = 20/c ^> 1. This free energy also exhibits 
a phase separation between a dense, highly connected 
phase in equilibrium with a dilute solution of beads. 

The mean field picture can be improved by using an 
exact mapping of the droplet-polymer system to a lattice 
gas model. Denoting the chemical potential of the beads 
as /if, and that of the polymers as the grand canonical 
partition function of the polymer-microemulsion system 
(with no loops) can be written as 



N p \ (P - 

^2Cl(N,P)e fibN (l-e^r P (5) 



9 = £"W*)E « 

N.P N„=0 p v ; 



N.P 



here £l(N, P) is the number of ways to arrange N drops 
so that P nearest neighbor pairs of drops exist in a 
given realization of the grand canonical ensemble {P is 
equal to \cN only in the mean field approximation). 
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Eq. JSJ) exactly corresponds to the grand canonical 
partition function of a lattice-gas with a Hamiltonian 
H = ln(l - e^) Y« l ,j > a i°~i + Mb Li Thus, the poly- 
mer chemical potential related to the number of the 
polymers in the system, plays the role of the interaction 
parameter that controls the phase separation. Note that 
this mapping can be extended off-lattice using similar, al- 
though more complicated, arguments. The very general 
nauture of this mapping pertains to any system of linked 
beads, and the results described in thsi Letter apply to 
many different experimental realizations, including DNA 
nanoparticles assemblies Q and antigen-antibodies im- 
munoessays Q 

Percolation: The transition from a fluid-like state to 
an elastic gel is related to the formation of an infinite 
network of droplets connected by polymers. This is de- 
scribed by site-bond percolation models that have been 
applied to study gelation ^(|. The vertices of a lattice 
are occupied with probability p s , and a bond can form be- 
tween two, occupied, nearest neighbors with probability 
Pb- Percolation is said to occur when an infinite cluster of 
occupied sites, connected by the bonds, is formed. In our 
system, the site-occupation probability p s can be identi- 
fied with the droplet concentration, c. The "bond occu- 
pation" probability can be identified as the probability 
that at least one polymer connects a given pair of beads. 
For a given numbers of chain ends per droplet, r, the av- 
erage number of bonding polymers per available nearest 
neighbor pair of drops is = r ( 1 ~ A ) where A is the av- 
erage fraction of looped polymers, Eq.Q. To a good ap- 
proximation, the probability p n that a given pair of near- 
est neighbor droplets is connected by n polymers, obeys 
the Poisson distribution, p n = ^-e~ rb . Thus, bond prob- 
ability is pi, = Y^=iPn = 1 — e~ Tb . Previous studies of 
percolation have shown |24j that the percolation line in 
the (pb,p s ) plane follows a power lawpj, = (p s )"pjj where 



a = — j^-p- and pi and p c s are the percolation thresholds 

for independently calculated bond and site percolation 
on the given lattice. Although lattice models for perco- 
lation cannot be applied in a quantitative manner to the 
continuum, polymer-droplet system they provide a qual- 
itative indication of the gelation transition as shown in 
Fig. 2. 

Simulations :To independently verify the predictions of 
the analytical model, we have performed the Monte 
Carlo simulations in the grand-canonical ensemble of 
both beads and polymers, with chemical potentials [ib 
and Hp, respectively. The simulated phase behavior also 
exhibits a first order transition between a dense con- 
nected network and a dilute phase. The detailed results 
of the simulations will be presented elsewhere pgj . 
The predicted phase diagram, summarized in Fig. 2, re- 
produces the experimentally observ ed p hase behavior of 
the model transient network of Ref.|2l|, shown in Fig. 1 
a). The spatial phase separation described in this Let- 
ter originates from a purely entropic effect, observed in 
other physical systems as well (23, |23| Namely, the loss 
of entropy due to spatial inhomogeneity and formation 
of the dense phase, is overcompensated by the increase 
in the polymer configurations. The simplicity of the ex- 
perimental realization which allows for easy control over 
the system parameters, the described system can be ex- 
tended to cover a wide range of equilibrium and non- 
equilibrium networks, in particular to examine the role 
of entropic effects in the formation of spatially inhomo- 
geneous structures in non-equilibrium networks |25l l26| 
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